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ABSTRACT 

We report precise Dopplcr measurements of seven subgiants from Keck Observatory. All seven 
stars show variability in their radial velocities consistent with planet-mass companions in Kcplerian 
orbits. The host stars have masses ranging from 1.1 < M+/Mq < 1.9, radii 3.4 < R±/Rq < 6.1, and 
metallicities —0.21 < [Fe/H] < +0.26. The planets are all more massive than Jupiter (Mp sin. i > 
1 Mj up ) and have semimajor axes a > 1 AU. We present millimagnitude photometry from the T3 
0.4 m APT at Fairborn observatory for five of the targets. Our monitoring shows these stars to 
be photometrically stable, further strengthening the interpretation of the observed radial velocity 
variability. The orbital characteristics of the planets thus far discovered around former A-type stars 
are very different from the properties of planets around dwarf stars of spectral type F, G and K, and 
suggests that the formation and migration of planets is a sensitive function of stellar mass. Three 
of the planetary systems show evidence of long-term, linear trends indicative of additional distant 
companions. These trends, together with the high planet masses and increased occurrence rate, 
indicate that A-type stars are very promising targets for direct imaging surveys. 
Subject headings: techniques: radial velocities — planetary systems: formation — stars: individual 
(HD 4313, HD 95089, HD 181342, HD 206610, HD 180902, HD 136418, HD 212771) 



1. INTRODUCTION 

The field of exoplanetary science recently reached a 
major milestone with the first direct-imaging detec- 
tions of planetary systems around main sequence stars0. 
iKalas et al.l (120081 ) detected a single planet-sized object 
with a semimajor axis a » 120 AU, orbiting just inside 
of the dust belt around the nearby, young A4V star Fo- 
malhaut. The young A5V dwarf star HR8799 is orbited 
by a system of three sub stellar objects with semimajor 
axes a = {24,38,68} AU (|Marois et al.ll200§) . These re- 
markable systems share a number of characteristics in 
common. Both host-stars are A-type dwarfs with stellar 
masses > 1.5 M Q surrounded by debris disks, the planets 
are super- Jupiters with masses < 3 Mj up , and the com- 
panions orbit far from their central stars with unexpect- 
edly large semimajor axes (ranging from 20 — 120 AU). 
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That both systems were discovered orbiting A stars 
might at first seem unlikely, given that A stars make 
up less than 3% of the stellar population in the Solar 
neighborhood and because the star-planet contrast ra- 
tios are unfavorable compared to systems with fainter, 
less massive central stars. However, in light of recent 
discoveries from Doppler-based planet searches of mas- 
sive stars it is becoming apparent that A dwarfs may 
in fact be ideal target stars for direct imaging surveys 
(lHatzes et al.l 12003 iSetiawan et al.l 120051: iReffert et al. 
2006 1 : [Sato et al.l 120071: iNiedzielski et al.l 120071: iLiu et al. 
20081 Zollinger et al.l 120091 ) . Measurements of the fre- 
quency of giant planets around the "retired" counter- 
parts of A-type dwarfs (subgiants and giants) have found 
that the occurrence of Jovian planets scales with stel- 
lar mass: A-type stars (M* > 1.5 M Q ) are at least 5 
times more likely than M dwarfs to harbor a giant planet 
(IJohnson et al.M2007at IBowler et al.ll2010l: IJohnson et al.l 
2009). And just like the current sample of imaged plan- 
ets, Doppler-de tected planets around retired A stars are 
more massive (jLovis fc Mayorl [2007) and orbit farther 
from their stars than do planet s found around Sun- 
like, F, G and K (FGK) dwarfs (IJohnson et al.l I2007bt 
ISato et al.ll200l . 

Indeed, there is strong evidence that the orbital char- 
acteristics of planets around A stars are drawn from 
a statistical parent population that is distinct from 
those of planets around FGK dwarfs. IBowler et al.1 
(2010) performed a statistical analysis of planets de- 
tected in the Lick Subgiants Survey, which com- 
prises 31 massive stars (M* > 1.5 M Q ) monitored 
for the past 5 years . The mass-period distribu- 
tion of exoplanets around FGK dwarfs is typically 
described by a double-power-law relationship, with 
the frequency of planets rising toward lower masses 
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and remaining flat in logarith mic semimajor-axis bins 
from ~ 0.05 AU to ~ 5 AU (iTabachnik fc Tremainel 
200 21 iLineweaver fc Gretherll2003t iCumming et al.ll200a 
Johnson! 12009). Based on the 7 planet detections from 
the Lick survey, Bowler et al. concluded that the power- 
law indices of the distribution of planets around A stars 
and Sun-like stars differ at the 4-cr level; the planets in 
their sample all have Mp sin i > 1.5 Mj up and none orbit 
within 1 AU. However, their small sample size precluded 
a determination of the exact shape of the mass-period 
distribution. Fortunately, given the 26±|% occurrence 
rate measured from the Lick survey, it will not take long 
to build a statistical ensemble comparable to the collec- 
tion of planets around less massive stars. 

To increase the collection of planets detected around 
massive stars, and to study the relationships among 
the characteristics of stars and the properties of their 
planets, we are conducting a survey of massive sub- 
giants at Keck and Lick Observatories. The de- 
creased rotation rates and cooler surface temperatures 
of these evolved stars make them much more ideal 
Doppler-survey targets compared to their m assive main- 
sequence prog enitors (|Galfand et al.l [2001 . The ob- 
served effects of stellar mass on the properties of plan- 
ets have i mportant implications for planet formation 
modeling (llda fc Linl 120051: Kennedy fc Kenvonl 120081: 
Kretke et al.ll2009t ICurridl2009t IDodson-Robinson et al l 
2009); the interpretation of observed st ructural fea- 
tures in the disks around massive sta rs ([Wvatt et al.l 
[1991 lQuillen|[200l iBrittain et aT] [20091 ): and the plan- 
ning of current and future planet sea rch efforts, such 
as th e Gemini Planet Imager (GPI; [Macintos h et al.l 
120081) . Spectro-Polari metric High-c o ntrast Exoplanets 
REsearch (SPHERE; IClaudi et al] [20061). the Near 
Infra red Coronographic Imager (NICI; IXrtigau et al.l 
|200"1 , and Project 1640 ([Hinklev et al.ll2008t ). Our Lick 
survey has resulted in the discovery of 7 new Jovian 
plan ets orbiting evolved stars more massive than the 
Sun (IJohnson et al.l l2006l I2007bl 12001 IPeek et al.lIlM 
IBowler et al.H2010D . In this contribution, we present the 
first seven planets discovered in the expanded Keck sur- 
vey. 

2. A DOPPLER SURVEY OF SUBGIANTS AT 
KECK OBSERVATORY 

2.1. Target Selection 

We are monitoring the radial velocities (RV) of a sam- 
ple of 500 evolved stars at Keck Observatory. The Keck 
program expands upon our Doppler survey of 120 sub- 
giants iatIjkk_^2bsCTyiyiO^ ongoing since 
2004 ([Johnson et al.lllMlPeek et al.|[2009t) . The stars in 
the Lick program have now been folded into the Keck tar- 
get list and that subset of brighter subgiants (V < 7.25) 
is currently monitored at both observatories. We began 
the Keck survey in 2007 April for the majority of our 
target stars, and a handful of stars were part of the orig- 
inal Keck planet se arch sample dating as far back as 1997 
([Marcv et al.l[2008l ). 

We selected the targets for the expanded Keck survey 
from the Hipparcos catalog based on the cri teria 1.8 < 
M v < 3.0, 8 < B - V < 1.1, and V < 8.5 ([ESAII1997I: 
van Lceuwcn 2007). We chose the red cutoff to avoid red 
giants, the majority of which are already monitored by 



other planet search progr ams and are known to exhibit 
veloc i ty jitter > 10 m s" 1 (iSato et al.ll2005l: iHekker et all 
2006; Ni edzielski et al.ll2009D . The lower M v restriction 
avoids Cepheid variables, and the upper limit excludes 
stars with masses less than 1.3 M when compared to 
th e Solar-metall i city ( [Fe/Hl — 0) stellar model tracks 
of iGirardi et al.l ((2002). We also excluded stars in the 
clump region (B - V > 0.8, M v < 2.0) to avoid the 
closely-spaced, and often overlapping mass tracks in that 
region of the theoretical H-R diagram. 

2.2. Stellar Properties 

Atmospheric parameters of the target stars are esti- 
mated from iodine-free, "template" spectra using the 
LTE spectro scopic analysis package S pectroscopy Made 
Easy (SME: IValenti fc Piskuno vlll996D. as descr i bed b y 
iValenti fc Fischerl (j2005| ) and iFischer fc Valentil (p005l ). 
To constrain the low surface gravit ies of the evolv e d star s 
we used the iterative scheme of IValenti et al.l (|2009l) . 
which ties the SME-derived value of log q to the gravity 
inferred from the Yonsei-Yale (Y 2 ; lYi et al.ll2004l) stellar 
model grids. The analysis yields a best-fit estimate of 
T c ff, log 3, [Fe/H], and V ro tsini. The properties of the 
majority of our targets from Lick and Keck are listed in 
the fourth edition of the Spectroscopic Properties of Cool 
Stars Catalog (SPOCS IV.; Johnson et al. 2010, in prep). 
We adopt the SME parameter uncertainties d escribed in 
the error analysis of IValenti fc Fischerl ([20051 ). 

The luminosity of each star is estimated from the 
appar ent V-band magnitude, the bolometric correc- 
tion dFlowerl 119961) . and the parallax from Hipparcos 
(|van Leeuwenl l2007f). From T e g and luminosity, we de- 
termine the stellar mass, radius, and an age estimate by 
associating those observed propertie s with a mode l from 
the Y 2 stellar interior calculations ([Yi et al.l 120041 ) . We 
also m easure the chromospheric emission i n the Ca II line 
cores (Wright et al .1120041 : llsaacsonl |2~009L providing an 
Shk value on the Mt. Wilson system , which we convert 
to \ogR' UK as per lNoves et alj if 19841 ). 

The stellar properties of the seven stars presented 
herein are summarized in Table [T] 

3. OBSERVATIONS AND ANALYSIS 
3.1. Keck Spectra and Doppler Analysis 

We obtained spectroscopic observations at Keck Obser- 
vatory using the HIRES spectrometer with a resolution 
of R k, 55, 000 with the B5 deck er (0"86 width) and red 
cross-disperser (|Vogt et alJll994f) . We use the HIRES ex- 
posure meter to ensure that all observations receive uni- 
form flux levels independent of atmospheric transparency 
variations, and to provide the photon- weighted exposure 
midpoint which is used for the barycentric correction. 
Under nominal atmospheric conditions, a V — 8 target 
requires an exposure time of 90 seconds and results in a 
signal-to-noise ratio of 190 at 5500 A. 

The spectroscopic observations are made through a 
temperature-controlled Pyrex cell containing gaseous io- 
dine, which is placed at the entrance slit of the spectrom- 
eter. The dense set of narrow molecular lines imprinted 
on each stellar spectrum from 5000 A to 6000 A provides 
a robust, simultaneous wavelength calibration for each 
observation, as well as information abou t the shape of the 
spectrometer's instrumental response ([Marcv fc Butlerl 
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1992). Doppler shifts are measured from e ach spectrum 
using the modeling procedure described bv lButler et al.l 
(1996). The instrumental uncertainty of each measure- 
ment is estimated based on the weighted standard devi- 
ation of the mean Doppler-shift measured from each of 
~ 700 2-A spectral chunks. In a few instances we made 
two or more successive observations of the same star and 
binned the velocities (in 2hr time intervals), thereby re- 
ducing the associated measurement uncertainty. 

3.2. Photometric Measurements 

We also acquired brightness measurements of five of 
the seven planetary candidate host stars with the T3 
0.4 m automatic photometric telescope (APT) at Fair- 
born Observatory. T3 observed each program star dif- 
ferentially with respect to two comparisons stars in 
the following sequence, termed a group observation: 
K,S,C,V,C,V,C,V,C,S,K, where if is a check (or sec- 
ondary comparison) star, C is the primary comparison 
star, V is the program (normally a variable) star, and 
S is a sky reading. Three V — C and two K — C dif- 
ferential magnitudes are computed from each sequence 
and averaged to create group means. Group mean dif- 
ferential magnitudes with internal standard deviations 
greater than 0.01 mag were rejected to eliminate the ob- 
servations taken under non-photometric conditions. The 
surviving group means were corrected for differential ex- 
tinction with nightly extinction coefficients, transformed 
to the Johnson system with yearly-mean transformation 
coefficients, and treated as single observations thereafter. 
The typical precision of a single group-mean observa- 
tion from T3, as measure d for pairs of con stant stars, 
is -0.004-0.005 mag (e.g. JHenrv eTaT]|205ol tables 2 & 
3). Further info rmation on the opera tion of the T3 APT 
can b e found in iHenrv et al.l (|1995bl lal) and lEaton et al.l 

(pomh . 

Our photometric observations are useful for eliminat- 
ing potential false positiv es fr om the sample of new plan- 
ets. lOueloz et all (|200lD and iPaulson et all (|2004l) have 
demonstrated how rotational modulation in the visibil- 
ity of starspots on active stars can result in periodic ra- 
dial velocity variations and potentially lead to erroneous 
planetary detections. Photometric results for the stars 
in the present sample are given in Table [5] Columns 
7-10 give the standard deviations of the V — C and 
K — C differential magnitudes in the B and V passbands 
with the 3cr outliers removed. With the exception of 
HD 206610, all of the standard deviations are small and 
approximately equal to the measurement precision of the 
telescope. 

For HD 206610, the standard deviations of the four 
data sets (V-C) B , (V-C) v , {K-C) B , and {K-C)y 
are all larger than 0.01 mag and indicate photometric 
variability. Periodogram analyses revealed that all four 
data sets have a photometric period of 0.09 day and an 
amplitude of — 0.03 mag. Thus, the variability must 
arise from HD 206610's primary comparsion star (C = 
HD 205318), which is included in all four data sets. Given 
its period, amplitude, and early-F spectral class, it is 
probably a new S Scuti star. We computed new differen- 
tial magnitudes for HD 206610 using the check star (K) 
to form the variable minus check data sets (V — K)b 
and (V — K)v- The standard deviations of these two 
data sets are 0.0066 and 0.0067 mag, respectively. 



Therefore, all five of the planetary candidate stars in 
Table [2] are photometrically constant to the approximate 
limit of the APT observations. The measured photomet- 
ric stability supports the planetary interpretation of the 
radial velocity variations. The two stars that we did not 
observe photometrically, HD 136418 and and HD 181342, 
have similar masses, colors, and activity levels (Tabled 
as the five stars we did observe and so are likely to be 
photometrically constant as well. 

3.3. Orbit Analysis 

For each star we performed a thorough search of the 
measured velocities for the best-fitting, single-planet Ke- 
plerian orbital model using the partial ly-linearized, least- 
square s fitting procedure described in iWright k, Howardl 
(2009) and implemented in the IDL package RVLItfl Be- 
fore searching for a best-fitting solution, we increased 
the measurement uncertainties by including an error con- 
tribution due to stellar "jitter." The jitter accounts for 
any unmodeled noise sources intrinsic to the star such 
as rotatio nal modulation o f surface inhomogeneities and 
pulsation feaar et al.lll998t I Wright! 120051 iMakarov et~aTI 
l2009HLagrange et alJl2010f >. 
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Fig. 1. — Top — Distribution of RVs for 72 standard stars, com- 
prising a total of 382 measurements. The dashed line shows the 
best-fitting Gaussian with a width a = 5.1 m s _1 . Bottom - 
Distribution of internal measurement uncertainties for 382 RV mea- 
surements. The median is 1.2 m s — 1 . Together with the distribu- 
tion of RVs in the top panel, this provides us with a jitter estimate 
of 5 m s" 1 , which we apply to all of the stars presented herein. 

We estimate the jitter for our subgiants based on the 
velocity variability of a sample of "stable" stars, for 

9 http://exoplanets.org/code/ 
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which we have obtained > 4 observations over a time 
span greater than 2 years. These stars do not show ev- 
idence of an orbital companion, except in a few cases 
where the stars exhibit a linear trend. In those cases 
we remove the trend using a linear fit and consider the 
scatter in the residuals. Figure [T] shows the distribution 
of RVs for all 72 stable stars, comprising 382 measure- 
ments. We fit a Gaussian function to the distribution 
with a width a = 5.1 m s _1 . The measurement uncer- 
tainties, shown in the lower panel, span 0.6-2.5 m s _1 , 
with a median value of 1.2 m s _1 . We subtracted this 
median internal error in quadrature from the measured 
width of the distribution of RVs to produce a jitter es- 
timate of 4.95 m s" 1 . In the analysis of each star's RV 
time series we round this value up and adopt a uniform 
jitter estimate of 5 m s _1 , which we add in quadrature to 
the measurement uncertainties before searching for the 
best-fitting orbit. 

After identifying the best-fitting model, we use a 
Markov-Chain Monte Carlo (MCMC) algor ithm to es- 
timate the p arameter uncertainties (See, e.g. lFordll2005l : 
iWinn et alll2007t ). MCMC is a Bayesian inference tech- 
nique that uses the data to explore the shape of the like- 
lihood function for each parameter of an input model. 
At each step, one parameter is selected at random and 
altered by drawing a random variate from a normal dis- 
tribution. If the resulting x 2 value for the new trial orbit 
is less than the previous x 2 value, then the trial orbital 
parameters are added to the chain. If not, then the prob- 
ability of adopting the new value is set by the ratio of the 
probabilities from the previous and current trial steps. If 
the current trial is rejected then the parameters from the 
previous step are adopted. 

We alter the standard deviations of the random pa- 
rameter variates so that the acceptance rates are be- 
tween 20% and 40%. The initial parameters are chosen 
from the best-fitting orbital solutions derived using the 
least-squares method described above, and each chain is 
run for 10 7 steps. The initial 10% of the chains are ex- 
cluded from the final estimation of parameter uncertain- 
ties to ensure uniform convergence. We verify that con- 
vergence is reached by running five shorter chains with 
10 6 steps and checking that the Gelman-Rubin statistic 
(|Gelman fc Ru bin 1992) for each parameter is near unity 
(< 1.02) and that a visual inspection of the history plots 
suggests stability. 

The resulting "chains" of parameters form the poste- 
rior probability distribution, from which we select the 
15.9 and 84.1 percentile levels in the cumulative distri- 
butions as the "one-sigma" confidence limits. In most 
cases the posterior probability distributions were approx- 
imately Gaussian. 

3.4. Testing RV Trends 

We use the Bayes i an In formation Criterion (BIC; 
ISchwarzl [19781 : ILiddlel 12001 and the MCMC posterior 
probability density functions (pdf) to deter mine whether 
there is evidence for a linear velocity trend ([Bowler et al.l 
120101) . The BIC rewards better- fitting models but penal- 
izes overly complex models, and is given by 

BIC = -21n£ max + fclnAr, (1) 
where £ max °c exp(— Xmin/ty is the maximum likelihood 



for a particular model with k free parameters and N 
data pointful. A difference of > 2 between BIC values 
with and without a trend indicates t hat there is s ufficient 
evidence for a more complex model (Kuha 2004). 

We also use the MCMC-derived pdf for the velocity 
trend parameter to estimate the probability that a trend 
is actually present in the data. If the 99.7 percentile 
of the pdf lies above or below m s" 1 yr _1 then we 
adopt the model with the trend. The BIC and MCMC 
methods yield consistent results for the planet candidates 
described in § [4] 

3.5. False-Alarm Evaluation 

For each planet candidate we consider the null- 
hypothesis that the apparent periodicity arose by chance 
from larger-than-expected radial velocity fluctuations 
and sparse sampling. We test this possibility by cal- 
culating the false-alarm probabili ty (FAP) based on 
the goodness-of-fi t statistic Ayf, ([Howard et al.l 120091 : 
iMarcv et ail 120051: ICummingj 12004ft , which is the differ- 
ence between two values of xt'- one from the single-planet 
Keplerian fit and one from the fit of a linear trend to the 
data. Each trial is constructed by keeping the times of 
observation fixed and scrambling the measurements, with 
replacement. We record the Axi value after each trial 
and repeat this process for 10,000 trial data sets. For 
the ensemble set we compare the resulting distribution 
of Ax 2 , to the value from the fit to the original data. 
The planets presented below all have FAP < 0.001, cor- 
responding to < 0.5 false alarms for our sample of 500 
stars. 

4. RESULTS 

We have detected seven new Jovian planets orbiting 
evolved, subgiant stars. The RV time series of each host- 
star is plotted in Figures [2HS1 where the error bars show 
the quadrature sum of the internal errors and the jit- 
ter estimate of 5 m s _1 , as described in § 13.31 The RV 
measurements for each star are listed in Tables [3H9[ to- 
gether with the Julian Date of observation and the in- 
ternal measurement uncertainties (without jitter). The 
best-fitting orbital parameters and physical character- 
istics of the planets are summarized in Table [101 along 
with their uncertainties. When appropriate we list notes 
for some of the individual planetary systems. 

HD 95089, HD 136418, HD 180902— The orbit models 
for these three stars include linear trends. We interpret 
the linear trend as a second orbital companion with a 
period longer than the time baseline of the observations. 

HD 181342— The time sampling of HD 181342 is 
sparser than most of the the other stars presented in this 
work. However, the large amplitude of the variations and 
observations clustered near the quadrature points result 
in a well-defined x 2 minimum in the orbital parameter 
space. The FAP for the orbit solution is 0.0064%. 

HD 212771 — This low-mass subgiant has a mass M* = 
1.15 Mq, indicating that it had a spectral type of early- 
G to late-F while on the main sequence. In addition 
to being one of our least massive targets, HD 212771 is 

10 The relationship between £ max and Xmin ' s om y valid under 
the assumption that the errors are described by a Gaussian, which 
is approximately valid for our analyses. 
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Fig. 2. — Relative RVs of HD4313 measured at Keck Obser- 
vatory. The error bars are the quadrature sum of the internal 
measurement uncertainties and 5 m s" 1 of jitter. The dashed line 
shows the best-fitting orbit solution of a single Keplerian orbit. 
The solution results in residuals with an rms scatter of 3.7 m s — 1 
and x/yj = 0.79. 



Fig. 4.— Relative RVs of HD 181342 measured at Keck Obser- 
vatory. The error bars are the quadrature sum of the internal 
measurement uncertainties and 5ms -1 of jitter. The dashed line 
shows the best-fitting orbit solution of a single Keplerian orbit. 
The solution results in residuals with an rms scatter of 4.7 m s — 1 
and w v? = 1.14. 
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Fig. 3.— Relative RVs of HD 95089 measured at Keck Obser- 
vatory. The error bars are the quadrature sum of the internal 
measurement uncertainties and 5 m s _1 of jitter. The dashed line 
shows the best-fitting orbit solution of a single Keplerian orbit plus 
a linear trend (dv/dt = 30.9±1.9 m s _1 yr — The solution results 
in residuals with an rms scatter of 5.7 m s _1 and \fx^ = 1.32. he 
lower panel shows the RVs with the linear trend removed. 



also one of the most metal-poor stars in the Keck sample 
[Fe/H] = -0.21. 

5. SUMMARY AND DISCUSSION 

We report the detection of seven new Jovian planets 
orbiting evolved stars. These detections come from the 
sample of subgiants that we are monitoring at Lick and 
Keck Observatories. The host-stars have masses in the 
range 1.15 M Q to 1.9 M , radii 3.4 < R±/R Q < 6.1, and 
metallicities -0.21 < [Fe/H] < +0.26. Five of the host- 
stars have masses M* > 1.5 M , and are therefore the 
evolved counterparts of the A- type stars. We also de- 
rived a jitter estimate for our sample of evolved stars 
and find that subgiants are typically stable to within 
5 m s" 1 . The observed jitter of subgiants makes them 
uniqu ely stable Doppler targets among mass ive, evolved 
st ars dFischer et alj|2003t iHekker et aLll2006ft . 

iBowler et al.l (|2010fl found that the minimum masses 
and semimajor axes of planets around A stars are very 
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Fig. 5.— Relative RVs of HD 206610 measured at Keck Obser- 
vatory. The error bars are the quadrature sum of the internal 
measurement uncertainties and 5ms -1 of jitter. The dashed line 
shows the best-fitting orbit solution of a single Keplerian orbit. 
The solution results in residuals with an rms scatter of 4.8 m s _1 
and \fxf. = 1.06. 



different from those of planets around FGK stars. Their 
findings suggest that the formation and migration mech- 
anisms of planets changes dramatically with increas- 
ing stellar mass. The planets reported in this work 
strengthen that conclusion. The five new planets we 
have discovered around stars with M± > 1.5 M Q all or- 
bit beyond 1 AU and have minimum masses Mp sin i > 
1 Mj up . These properties contrast with those of plan- 
ets orbiting less massive stars, which have a nearly flat 
distribution in log a f rom a = 0.05 AU to a = 1 AU 
(jCumming et al.|[2008l) . and a steeply rising mass func- 
tion with din N/d In Mp = —1.4. Successful theories of 
the origin and orbital evolution of giant planets will need 
to account for the discontinuity between the distribu- 
tions of orbital p arameters for planets around Sun-like 
and A-type stars (Kennedy fc Kenyonl2008t lCurriell2009t 
IKretke et al.ll2009D . 

The abundance of super- Jupiters (Mp s'mi > 1 Mj up ) 
detected around massive stars bodes well for future 
direct-imaging surveys. In addition to harboring the 
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Fig. 6.— Relative RVs of HD 180902 measured at Keck Obser- 
vatory. The error bars are the quadrature sum of the internal 
measurement uncertainties and 5 m s" 1 of jitter. The dashed line 
shows the best-fitting orbit solution of a single Keplerian orbit plus 
a linear trend (dv/dt = 135.4 ± 3.5 m s — 1 yr _1 ). The solution re- 
sults in residuals with an rms scatter of 3.3 m s — 1 and \J x% = 0.98. 
The lower panel shows the RVs with the linear trend removed. 
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Fig. 7.— Relative RVs of HD 136418 measured at Keck Obser- 
vatory. The error bars are the quadrature sum of the internal 
measurement uncertainties and 5 m s _1 of jitter. The dashed line 
shows the best-fitting orbit solution of a single Keplerian orbit plus 
a linear trend (dv/dt = —5.3 ± 1.7 m s _1 yr — 1 ). The solution re- 
sults in residuals with an rms scatter of 5.0 m s — 1 and y/xl = 1.15. 
HD 136418 has a mass M* = 1.33 Mq, making it a former F-type 
star. 

massive planets that are predicted to be the most easily 
detectable in high-contrast images, A-type dwarfs have 
the added benefit of being naturally young. A 2 M© 
star has a main-sequence lifetime of only ~ 1 Gyr, which 
means that Jovian planets in wide orbits will be young 
and thermally bright. Three of the planets in Table [TU] 
show linear velocity trends indicative of additional long- 
period companions. These linear trends provide clear 



markers of massive objects in wide orbits around nearby 
stars, and therefore warrant additional scrutiny from RV 
monitoring and high-contrast imaging. 



We thank the many observers who contributed to the 
observations reported here. We gratefully acknowledge 
the efforts and dedication of the Keck Observatory staff, 
especially Grant Hill, Scott Dahm and Hien Tran for 
their support of HIRES and Greg Wirth for support 



100 



50 



-50 



-100 



" P = 373.3 d 


i ... | 

HD 212771 


K = 58.2 m s _1 


Keck Obs. 


.8 = 0.11 






/ \ 


f 




- V \ 


/ W - 




+ 


MpSini = 2.3 M Jup 






1 ... 1 



2007 



2008 



2009 



2010 



Year 



Fig. 8. — Relative radial velocity measurements of HD 212771. 
The error bars are the quadrature sum of the internal measure- 
ment uncertainties and 5 m s _1 of jitter. The dashed line shows 
the best-fitting orbit solution of a single Keplerian orbit. The so- 
lution results in residuals with an rms scatter of 5.8 m s _1 and 
^/xl = 1-29. HD 212771 has a mass M* = 1.15 M , indicating 
that it was either an early-G or late-F star when it was on the main 
sequence. 
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TABLE 1 
Stellar Parameters 



— 

Parameter 


HD 4313 


HL) 95089 


HD 181342 


HI) 2U661U 


HD 18U9U2 


HD 136418 


ML) 21277 1 


V 


7.83 


7.92 


7.55 


8.34 


7.78 


7.88 


7.60 


B — V 


0.96 


0.94 


1.02 


1.01 


0.94 


0.93 


0.88 


Distance (pc) 


137 (14) 


139 (16) 


110.6 (7.5) 


194 (36) 


110 (10) 


98.2 (5.6) 


131 (14) 


M v 


2.2 (0.3) 


2.1 (0.3) 


2.2 (0.2) 


2.2 (0.4) 


2.5 (0.3) 


2.7 (0.2) 


2.2 (0.3) 


[Fe/H] 


+0.14 (0.03) 


+0.05 (0.03) 


+0.26 (0.03) 


+0.14 (0.03) 


0.04 (0.03) 


-0.07 (0.03) 


-0.21 (0.03) 


T eff (K) 


5035 (44) 


5002 (44) 


5014 (44) 


4874 (44) 


5030 (44) 


5071 (44) 


5121 (44) 


Vrot sm i (km s 1 ) 


2.76 (0.5) 


2.74 (0.5) 


3.04 (0.5) 


2.57 (0.5) 


2.88 (0.5) 


0.17 (0.5) 


2.63 (0.5) 




3.4 (0.06) 


3.4 (0.06) 


3.4 (0.06) 


3.3 (0.06) 


3.5 (0.06) 


3.6 (0.06) 


3.5 (0.06) 


M, (M ) 


1.72 (0.12) 


1.58 (0.11) 


1.84 (0.13) 


1.56 (0.11) 


1.52 (0.11) 


1.33 (0.09) 


1.15 (0.08) 


R* (R©) 


4.9 (0.1) 


4.9 (0.1) 


4.6 (0.1) 


6.1 (0.1) 


4.1 (0.1) 


3.4 (0.1) 


5.0 (0.1) 


L» (R Q ) 


14.1 (0.5) 


13.5 (0.5) 


12.0 (0.5) 


18.9 (0.6) 


9.4 (0.5) 


6.8 (0.5) 


15.4 (0.5) 


Age (Gyr) 


2.0 (0.5) 


2.5 (0.9) 


1.8 (0.4) 


3(1) 


2.8 (0.7) 


4(1) 


6(2) 


Shk 


0.12 


0.13 


0.12 


0.14 


0.15 


0.14 


0.16 




-5.27 


-5.22 


-5.31 


-5.23 


-5.14 


-5.19 


-5.09 



TABLE 2 

Summary of Photometric Observations From Fairborn Observatory 



Program 


Comparison 


Chock 


Date Range 


Duration 




"(V - C) B 


a(V - C)y 


a(K - C) B 


<r(if - C)y 




Star 


Star 


Star 


(HJD - 2,440,000) 


(days) 


™obs 


(mag) 


(mag) 


(mag) 


(mag) 


Variability 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


O) 


(10) 


(11) 


HD 4313 


HD 4627 


HD 4526 


54756-55222 


466 


188 


0.0039 


0.0046 


0.0038 


0.0056 


Constant 


HD 95089 


HD 94401 


HD 93102 


55139-55237 


98 


58 


0.0056 


0.0059 


0.0045 


0.0064 


Constant 


HD 180902 


HD 179949 


HD 181240 


55122-55126 


4 


3 


0.0032 


0.0031 


0.0010 


0.0056 


Constant? 


HD 206610 


HD 205318 


HD 208703 


54756-55170 


414 


88 


0.0183 a 


0.0140 3 " 


0.0172 


0.0124 


Constant 


HD 212771 


HD 212270 


HD 213198 


55119-55170 


51 


60 


0.0043 


0.0048 


0.0051 


0.0065 


Constant 



a Comparison star HD 205318 is variable in brightness, so wc recomputed the standard deviations of the program star from the (V — K)b 
and (V — K)v differential magnitudes and find them to be 0.0066 and 0.0067 mag in B and V, respectively. 



TABLE 3 
Radial Velocities for HD 4313 



JD 


RV 


Uncertainty 


-2440000 


Cm s- 1 ) 


(ms-1) 


14339.932 


23.92 


1.57 


14399.842 


21.61 


1.59 


14456.806 


-30.29 


1.61 


14675.006 


8.82 


1.71 


14689.004 


13.91 


1.60 


14717.945 


24.71 


1.53 


14722.895 


31.39 


1.60 


14779.854 


8.37 


1.70 


14790.889 


1.82 


1.64 


14805.807 


-18.24 


1.54 


14838.768 


-46.70 


1.61 


14846.745 


-54.50 


1.69 


14867.754 


-61.16 


3.56 


14987.118 


-25.11 


1.69 


15015.049 


-3.83 


1.54 


15016.081 


-1.47 


1.41 


15027.089 


8.26 


1.61 


15049.038 


22.34 


1.56 


15076.091 


22.90 


1.60 


15081.091 


21.97 


1.55 


15084.143 


15.79 


1.60 


15109.955 


15.55 


1.63 


15133.975 


0.00 


1.60 


15169.860 


-23.21 


1.57 


15187.855 


-42.72 


1.55 


15198.771 


-51.25 


1.55 


15229.722 


-63.75 


1-43 


15250.713 


-67.93 


1.50 



TABLE 4 
Radial Velocities for HD 95089 



JD 


RV 


Uncertainty 


-2440000 


(m a" 1 ) 


(m a" 1 ) 


14216.851 


-11.35 


1.26 


14429.126 


3.46 


1.35 


14543.962 


27.54 


1.44 


14602.807 


-1.61 


1.19 


14603.793 


7.19 


1.39 


14808.056 


-21.62 


1.25 


14839.102 


-20.36 


1.49 


14847.054 


-10.63 


1.52 


14865.063 


-12.19 


1.27 


14929.816 


4.82 


1.47 


14956.912 


21.43 


1.36 


14983.770 


14.90 


1.40 


14984.829 


17.03 


1.50 


14987.835 


29.57 


1.28 


15015.761 


30.40 


1.33 


15135.149 


-0.70 


1.42 


15164.109 


0.53 


1.34 


15172.146 


0.00 


1.45 


15173.095 


-4.48 


1.47 


15188.103 


-2.29 


1.42 


15232.135 


-20.55 


1.49 


15252.034 


-30.59 


1.53 



TABLE 5 
Radial Velocities for HD 181342 



JD 


RV 


Uncertainty 


-2440000 


(m s-1) 


(ms-1) 


14339.768 


2.69 


1.28 


14399.741 


0.00 


1.33 


14634.063 


-77.87 


1.43 


14641.000 


-81.00 


1.44 


14674.973 


-91.85 


1.32 


14808.687 


-66.72 


1.50 


14957.027 


8.89 


1.28 


14964.119 


10.11 


1.29 


14984.083 


21.20 


1.31 


14985.112 


21.37 


1.35 


15015.014 


9.33 


1.22 


15016.963 


11.24 


1.41 


15026.967 


7.63 


1.38 


15042.963 


-1.18 


1.38 


15109.749 


-27.32 


1.34 


15135.743 


-22.92 


1.38 


15169.686 


-43.74 


1.42 



TABLE 6 
Radial Velocities for HD 206610 



,JD RV Uncertainty 



-2440000 


(m s-1) 




(m 


14313 


980 


22.85 


1 


.48 


14339 


843 


2.84 


1 


.25 


14399 


759 


0.00 


1 


48 


14674 


964 


39.03 


1 


12 


14717 


919 


54.16 


1 


20 


14779 


822 


59.68 


1 


.36 


14790 


746 


52.55 


1 


42 


14807 


783 


55.93 


1 


.57 


14956 


102 


11.16 


1 


28 


14964 


118 


6.98 


1 


.26 


14984 


082 


-5.90 


1 


27 


14985 


111 


0.42 


1 


.43 


14986 


113 


-1.17 


1 


.27 


15015 


022 


-11.00 


1 


.36 


15015 


956 


-8.62 


1 


.37 


15019 


066 


-1.59 


1 


.47 


15029 


011 


3.61 


1 


.46 


15043 


058 


-10.34 


1 


31 


15077 


058 


-13.95 


1 


.40 


15083 


055 


-13.70 


1 


11 


15135 


761 


-23.43 


1 


13 


15169 


700 


-25.24 


1 


.46 


15187 


696 


-7.58 


1 


50 


15188 


691 


-15.51 


1 


.40 
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TABLE 7 
Radial Velocities for HD 180902 



JD 


RV 


Uncertainty 


-2440000 


Cm s" 1 ) 


(ms-1) 


14339.767 


40.51 


1.17 


14399.740 


19.80 


1.37 


14634.062 


-4.44 


1.37 


14640.998 


-14.25 


1.59 


14674.972 


5.16 


1.51 


14808.686 


37.39 


1.56 


14957.026 


0.00 


1.44 


15015.013 


-15.62 


1.39 


15026.965 


-23.67 


1.33 


15109.748 


-11.60 


1.45 


15135.741 


-1.92 


1.33 


15169.685 


14.50 


1.51 



TABLE 8 
Radial Velocities for HD 136418 



JD 


RV 


Uncertainty 


-2440000 


(m b-1) 


(m s" 1 ) 


14216.954 


25.79 


1-27 


14345.803 


27.56 


1-50 


14549.041 


-23.63 


1-97 


14637.956 


0.03 


1.26 


14674.797 


0.00 


1.32 


14839.165 


0.97 


1.21 


14927.948 


-40.01 


1.66 


14955.930 


-39.28 


1.38 


15014.866 


-43.17 


1-21 


15016.982 


-37.15 


1-13 


15028.959 


-31.44 


0.96 


15041.838 


-31.07 


1-26 


15042.875 


-32.70 


1.33 


15082.729 


-28.68 


1.46 


15106.719 


-13.60 


1.43 


15111.702 


-2.58 


1.49 


15188.168 


44.27 


1.27 


15189.147 


31.37 


1.37 


15192.148 


43.12 


1.31 


15197.173 


41.88 


0.99 


15231.100 


38.97 


1.41 


15232.083 


40.97 


1.32 


15252.046 


30.85 


1.26 


15256.997 


34.86 


1.34 



TABLE 9 
Radial Velocities for HD 212771 



JD 


RV 


Uncertainty 


-2440000 


(m b-1) 


(ms-l) 


14339.830 


-20.86 


2.03 


14399.773 


9.64 


1.79 


14456.785 


45.09 


1.68 


14675.025 


-13.62 


1.77 


14807.780 


38.82 


1.84 


14964.124 


49.23 


2.70 


14984.088 


34.58 


2.51 


14985.116 


20.73 


2.62 


14986.117 


29.29 


2.77 


15015.023 


-10.59 


2.63 


15016.078 


-13.55 


1.63 


15019.080 


-6.43 


1.88 


15027.008 


-8.30 


1.80 


15043.060 


-17.15 


1.96 


15076.073 


-34.91 


1.91 


15083.064 


-21.87 


1.74 


15134.930 


0.00 


1.74 


15135.765 


-8.56 


1.52 


15169.702 


16.03 


1.77 


15188.695 


44.57 


1.77 
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TABLE 10 
Orbital Parameters 



Parameter 


HD 4313 b 


HD 95089 b 


HD 181342 b 


HD 206610 b 


HD 180902 b 


HD 136418 b 


HD 212771b 


Period (d) 


356.0 (2.6) 


507 (16) 


663 (29) 


610 (13) 


479 (13) 


464.3 (3.2) 


373.3 (3.4) 


T p a (JD) 


2454804 (80) 


2454983 (90) 


2454881 (50) 


2454677 (20) 


2454690 (130) 


2455228 (10) 


2454947 (40) 


Eccentricity 


0.041 (0.037) 


0.157 (0.086) 


0.177 (0.057) 


0.229 (0.058) 


0.09 (0.11) 


0.255 (0.041) 


0.111 (0.060) 


K (m s" 1 ) 


46.9 (1.9) 


23.5 (1.9) 


52.3 (3.7) 


40.7 (1.9) 


30.7 (3.7) 


44.7 (1.9) 


58.2 (7.8) 


uj (deg) 


76 (80) 


317 (60) 


281 (30) 


296 (10) 


300 (100) 


12 (10) 


55 (40) 


Mpsini (M Jup ) 


2.3 (0.2) 


1.2 (0.1) 


3.3 (0.2) 


2.2 (0.1) 


1.6 (0.2) 


2.0 (0.1) 


2.3 (0.4) 


a (AU) 


1.19 (0.03) 


1.51 (0.05) 


1.78 (0.07) 


1.68 (0.05) 


1.39 (0.04) 


1.32 (0.03) 


1.22 (0.03) 


Linear trend (m s _1 yr — 1 ) 


(fixed) 


30.9 (1.9) 


(fixed) 


(fixed) 


135.4 (3.5) 


-5.3 (1.7) 


(fixed) 


rms (m s~ ) 


3.7 


5.7 


4.7 


4.8 


3.3 


5.0 


5.8 


Jitter (m s _1 ) 


5.0 


5.0 


5.0 


5.0 


5.0 


5.0 


5.0 


n/x? 


0.79 


1.32 


1.14 


1.06 


0.98 


1.15 


1.29 


N obs 


29 


22 


17 


24 


12 


25 


20 



a Time of pcriastron passage. 

b In cases where the eccentricity is consistent with e — 0, we quote the 2er upper limit from the MCMC analysis. 



